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Continuum and emission-line variability of active galactic nuclei provides a power- 
ful probe of microarcsecond scale structures in the central regions of these sources. 
In this contribution, we review basic concepts and methodologies used in analyz- 
ing AGN variability. We develop from first principles the basics of reverberation 
mapping, and pay special attention to emission-line transfer functions. We discuss 
application of cross-correlation analysis to AGN light curves. Finally, we provide 
a short review of recent important results in the field. 



1 Introduction 

The study of multiwavelength variability of active galactic nuclei (AGNs) is 
now a major branch of the field, enabled largely by the availability of suitable 
facilities for long-term studies of faint sources at many wavelengths; some of 
the scientific arguments for development of these facilities have been based on 
exactly these programs, which are now leading to improved understanding of 
the AGN phenomenon. It is only within the last few years that the evidence for 
supermassive black holes in both active and non-active galaxies has gone from 
circumstantial to compelling, and the potentially most powerful technique for 
measuring black-hole masses in AGNs is through study of broad emission- 
line variability. The existence of accretion disks in AGNs is still far from 
proven, but the evidence for them is improving, again as a result of variability 
studies. Multiwavelength monitoring observations are beginning to show the 
relationships between variability in different bands, and the hope is that once 
the phenomenology is better known, our understanding of the physics will 
follow. 

In this contribution, we will cover the basic characteristics of AGN vari- 
ability and provide what we hope is some relevant historical background. 
Because of the current importance of emission-line variability studies, we will 
develop the theory of reverberation mapping from first principles. One of 
the most powerful and widely used tools in the analysis of emission-line and 
continuum variability data is the technique of cross-correlation, and we will 
therefore describe in some detail the application of this method to AGN data. 
Throughout this contribution, we will concentrate almost exclusively on non- 
blazar AGNs, those for which we believe most of the observed UV/optical 
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emission originates in an accretion disk rather than in a relativistic jet: many 
of the techniques described here are also applicable to blazars, however. For 
a fairly recent discussion of blazar variability results, we refer the reader to 
the fine review by Ulrich et al.E£l. We emphasize that we intend for this con- 
tribution to be primarily instructional; this should not be misconstrued as a 
comprehensive review of the state of the field. Our intent is to provide both 
students and researchers who already have some familiarity with AGNs with 
enough background to read critically the current literature on AGN variability 
and understand the strengths and weaknesses of the method. 

2 Background and Basic Phenomenology 

AGNs show flux variations over the entire electromagnetic spectriixcL-. Indeed, 
variability was one of the first recognized properties of quasarsEil'Eii Early 
investigations established that significant variations ( £ 0.1 mag) in the optical 
brightness of quasars could occur on time scales as short as days. 

Detection of rapid variability in quasars was a remarkable discovery at the 
time because it implies that the size of the continuum-emitting region must 
be of order light days (1 lt-day = 2.6 x 10 15 cm), based on source coherence 
arguments: for a source to vary coherently, the entire emitting region must 
be causally connected, which implies a maximum size for the source based 
on light-travel time. Suppose, for example, that the brightness of the source 
doubles in a week; we can immediately conclude that the emitting region must 
be no larger than one light week in radius on the basis of causality. We might 
suppose that in fact there are multiple emitting regions varying at random, 
but the number of such regions must be limited or the stochastic variations 
would be averaged out. Individually, then, the purported independent regions 
face a similar size limit from causality, and the conclusion that the emitting 
regions are very small cannot be avoided. Historically, this is the quasar 
problem: how can so much energy, the equivalent of as much as trillions of 
stars, be produced in a region that is about the size of the Solar System? 

The detection of quasar variability was a critical part of the argument 
that AGNs are powered by supermassive black holes. The original arguments 
for supermassive black holes in AGNs were based on mass constraints from 
the Eddington limit and size constraints from variability!^! The Eddington 
limit is the requirement that gravitational forces on an ionized gas exceed 
outward radiation pressure, which translates to a requirement that 

M > M Edd = 8 x 10 5 ( -) M . (1) 

\ 10 44 ergs s 1 / 
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Rapid variations, in some cases on time scales as short as a day, require an 
emitting region less than a light day in radius, which for a 10 44 ergs s _1 AGN 
corresponds to ~ 10 4 i? grav , where i? gra v = GM/c is the gravitational radius. 
Later arguments for supermassive black-holes in AGNs focussed on how to 
power AGNs by gravitational accretionEil. The deep gravitational potential 
leads to an accretion disk that radiates most strongly across the UV/optical 
spectrum, and for AGN masses above the Eddington limit, the thermal emis- 
sion should peak in the near UV. Indeed, then, the broad UV/optical feature 
known as the "big blue bump" can plausibly be identified with accretion-disk 
emission. Furthermore, the intense magnetic fields expected in disks could 
provide a mechanism for jet collimation. 



2.1 Basic Characteristics of Variability 

AGNs have been found to be variable at all wavelengths at which they have 
been observed. The variations appear to be aperiodic and have variable ampli- 
tude. While variability in high-luminosity AGNs (quasars) was reported soon 
after their discovery, variability in lower-luminosity AGNs (Seyfert galaxies) 
was not reportedEj until 1967, and was less dramatic. The reason for this 
is probably quite simple: most of the quasars that were monitored are now 
known to be the jet-dominated sources known as "blazars". i.e., BL Lac ob- 
jects and optically violent variables (OVVs). The optical identifications of 
quasars were based on coincidence with radio-source positions, which natu- 
rally led to biases towards radio-loud quasars, blazars in particular. Indeed, 
the original arguments about size and variability time scales in retrospect 
apply to the jets, not necessarily what we now identify as thermally emit- 
ting accretion disks. Nevertheless, the original conclusions about AGN sizes 
proved to be generally correct for both blazars and non-blazars. 

UV / Optical Variability. Fig. ^ shows a light curve for a typical Seyfert 1 
galaxy, NGC 5548, which will serve as a continuing example through this 
chapter as it is one of the best-studied objects of this class. While no periodic 
behavior has been identified, there are some basic parameterizations that allow 
us to characterize the variability. An example is shown in Fig. ^| We compare 
each flux measurement at an arbitrary time U with the flux at every later time 
tj . The quantity we show is the flux ratio as a function of the time interval 
between observations At = tj — U, i.e., Alogi 7, = log F(tj)/F(ti), where in 
each case the contaminating flux due to starlight in the host galaxy (as shown 
in Fig. |j]) has been first removed. Fig. § shows that NGC 5548 shows little 
variability on time scales shorter than a few days, but on time scales of several 
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Figure 1. Optical (5100 A) light curve of NGC 5548 from late 1988 to late 1996. The 
horizontal line indicates the constant contribution from starlight through the standard 
aperture used, which projects to 5" X 7".5. From Peterson et al. 71 ©1999 AAS. 



weeks or months, very large variations can be observed. We note in passing 
that the quantity shown here is closely related to the "structure function". 
The structure function is simply the mean absolute value of 2. 5 A log F(At), 
i.e., the mean difference in magnitudes between observations separated by 
times At. 

A common parameter to characterize variability is the mean fractional 
variation, 



F = 



Vct 2 - g 



where the quantities are (a) the mean flux for all N observations, 



N 



</> = £/* . 



(2) 



(3) 
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Figure 2. Difference in opticalflux as a function of interval between observations At for 
the light curve shown in Fig. nj. The starlight contribution as shown in Fig. hi has been 
subtracted off. 



(b) the variance of the flux (as observed) , 
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and (c) the mean square uncertainty of the fluxes 




N 



(5) 



Even if the continuum is constant, there will still be apparent flux variations 
simply due to measurement uncertainties (noise) ; the virtue of F var is that it 
adjusts the fractional variation downwards to account for the effect of ran- 
dom errors. The parameter F vax is thus sometimes referred to as the "excess 
variance." In Fig. [j], we show F val as a function of time interval UV and 
optical variations in NGC 5548. This shows that UV continuum variations 
are typically around 10-20% on time scales of about a month. The fractional 
variations in the optical are less pronounced, at least in part because of con- 
tamination of the fluxes by a large constant contribution from starlight in the 
host galaxy. 

X-Ray Variability. Rapid _£-ray variability is a hallmark of AGNs (for a re- 
view, see Mushotzky et al.Ej). It is of historical importance since it effectively 
eliminates alternatives to supermassive black holes, e.g., massive stars or star- 
bursts, as competing explanations for the high luminosities of AGNs. X-rays 
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Figure 3. Variability parameter -F var from ultraviolet and optical continuum measurements 
of NGC 5548, as a function of interval between observations; the optical data used are those 
shown in Fig. hi Courtesy of S. Collier. 



are expected to arise near the event horizon of the black hole, so the shortest 
time-scale X-ray variability is expected on a few times the crossing time 

Crossing = = 10 ( 1Q6 S. (6) 

Prior to the 1990s and the advent of RXTE, the best X-ray monitoring 
data was from EXOSAT, which was in a high-Earth orbit that allowed up to 80 
hours of uninterrupted observations. The EXOSAT "long-looks" at variable 
AGNs established that variability is most rapid in low-luminosity systems. 

A useful way to characterize variability is in terms of the "power density 
spectrum" (PDS), which is the product of the Fourier transform of the light 
curve and its complex conjugate. The PDS for AGNs is often parameterized 
as a power law, 

p(f) = r a ■ (7) 

The EXOSAT dataS showed that AGN X-ray variations can be characterized 
by PDS indices in the range 1 a 2 over time scales of hours to months. 
The total power in the variations is given by integrating the PDS over all 
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frequencies. Thus, the PDS must turn over at low frequencies (i.e., a as 
defined above must become less than unity) to prevent divergence in the total 
power. Such breaks in the PDS are observed in stellar-mass X-ray sources, and 
the turnover frequency correlates inversely with mass, though the fundamental 
reason for this is not understood. The basic idea is that the mass of the black 
hole can be inferred, since Mbh oc i? gra v oc Crossing- If we scale AGNs relative 
to stellar-mass systems (which have turnover frequencies / w 0.1 Hz), we 
expect that the turnover frequencies for AGNs will occur at about 10~ 8 Hz. 
In only one case, NGC 3516, has thep. been a plausible detection of the 
turnover frequency. Edelson & NandraEil find that for this Seyfert 1 galaxy 
the turnover frequency is / w 4 x 10~ 7 Hz, which corresponds to a time scale 
of about one month. The mass inferred, again scaling relative to stellar-mass 
systems, is in the range 1O 6 -1O 7 M0. 

Periodicities in X-ray light curves, which might reflect orbital or preces- 
sion periods, have been searched for, but never found. As a historical-footnote, 
however, it is worth mentioning that one such detection was claimecO, namely 
a 12,000 s period in NGC 6814. However, ROSAT observations revealed that 
the variable source is in fact a foreground Galactic binary in the same field as 
the AGNB 

2.2 Origin of the Variations 

At a fundamental level, the physical origin of variations is not known, although 
accifition-disk instabilities are probably involved. For example, Kawaguchi et 
al.E3 show that continuum variations with a PDS of the form of Eq. (Q) can 
be explained by magnetohydrodynamic instabilities, specifically disconnection 
events, within the disk. Variable accretion rates have also been considered. 
In some specific atypical cases, variations have been attributed to variable 
obscuration of the nuclear source and to microlensing due to stars in the host 
galaxy. 

There are a number of important physical time scales that might be asso- 
ciated with variability We include them here in the convenient form given by 
Edelson & NandracJ. The first of these is the crossing time already mentioned 
m Eq. (§), which we rewrite here as 



which is the time it takes a radiative signal to cross the X-ray emitting region 
(assumed to be at r ^ 10i? grav ). Here M7 is the black-hole mass in units of 
10 7 M Q . Variations might also be expected on the time scale of the orbital 
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period, 



^orbital — 0.33M7 




(9) 



Thermal instabilities might also cause variations on the time scale for their 
development, 



where a is the viscosity parameter. Mechanical instabilities may propagate as 
acoustic waves, which will travel at the sound speed and thus cross the disk 
on a time scale 



where H is the disk thickness. And finally, the time scale over which the 
effects of variations in the accretion rate will propagate through the disk is 
given by the drift speed, 



2.3 Blazar Variability 

Variability properties of blazars are distinct from those of other non-beamed 
AGNs and should be mentioned separately. Blazars are characterized by ex- 
treme variability at all wavelengths. Unlike radio-quiet quasars or Seyfert 
galaxies, significant infrared and radio continuum variability is observed in 
blazars. Furthermore, the polarization of the continuum radiation is also sig- 
nificant (i.e., greater than a few percent), and the degree of polarization and 
amplitude of variability correlate with luminosity, which is the opposite case 
for non-blazars. Blazars are also the only sources detected at TeV energies, 
and the TeV fluxes can vary by as much as a factor of 10 in one day. All of 
these properties indicate that the continuum is dominated by emission from 
relativistic jets, as these characteristics suggest a non-thermal (synchrotron 
or inverse Compton) origin. 

2.4 Emission- Line Variability 

The broad emission lines in AGN spectra can vary both in flux and in profile. 
Over time scales of months and years, the changes can be very dramatic, but 
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on shorter time scales they are more subtle^. The first detection of emission- 
line variations was by Andrillat & Souflrincl, based on photographic snectca 
of the Seyfert 1 galaxy NGC 3516. There were a few subsequent reportsEJ'El, 
but these cases seemed to be widely regarded as "curiosities" that did not 
generate much follow-up work. The basic problem was that only very large 
changes could be detected photographically or with the intensified television- 
type scanners that were commonly used in AGN spectroscopy from the mid- 
1970s to mid-1980s. In the few cases where clear variations were detected, the 
changes were often dramatic; there were sometimes claims of Seyferts changing 
"type" as broad components of emission lines appeared or disappeared. 

Pronounced variability of broad emission-line profiles was detected in the 
early 1980s by a number of investigators. Profile variations were originally 
thought to be due to excitation inhomogeneities; excitation pulses propagat- 
ing through a broad-line region (BLR) with an ordered velocity field would 
produce features that could propagate across the profile with time. This con- 
cept led to the development of reverberation mapping, which is described in 
detail in Sec. 0. PetersonEil reviews early work on emission-line variability. 

A useful way to isolate the variable part of an emission line is shown in 
Fig. |. The upper panel shows a mean spectrum formed from 34 individual 
HST spectra of NGC 5548. The lower panel shows the "root-mean square" 
(rms) spectrum which is formed from the same data simply by computing the 
rms flux at each wavelength. Constant features, such as narrow emission lines 
and host-galaxy flux, do not appear in the rms spectrum. 

2.5 The First Monitoring Programs 

The early 1980s saw the first attempts to monitor the UV/optical continuum 
and emission-line variations in Seyfert 1 galaxies. There were two reasons this 
happened when it did. First, there was a realization that variability afforded 
a powerful probe of the structure and kinematics of AGNs on projected scales 
of microarcseconds. Variability was recognized as an important new tool 
with which to study enigmatic quasars. Second, the right technology for 
such investigations became widely available: it was possible to attempt such 
programs on account of (a) IUE, which for the first time allowed precision UV 
spectroscopy of low-redshift extragalactic objects, and (b) the proliferation of 
linear electronic detectors (first Reticons and Image Dissector Scanners, and 
later CCDs) on moderate-size (l-2m) ground-based telescopes. 

One of the first significant monitoring programs was a multiple-year IUE- 
based program on NGC 4151, which wasxarried out by a European consortium 
led by M.V. Penston and M.-H. UlricbEa. Ultraviolet spectra were obtained 
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Figure 4. The top panel shows the mean spectrum computed from 34 HST spectra of the 
variable Seyfert 1 galaxy NGC 5548 40 . The lower panel shows the rms spectrum based on 
variations around this mean. The rms spectrum thus isolates the variable components of 
the spectrum. Fluxes are in units of 10~ 15 ergs s — 1 cm~ 2 A -1 . 



with a typical sampling interval of 2-3 months. The program showed that 
variations in the UV and optical continua were closely coupled. It also revealed 
that the emission- line flux variations are correlated with continuum variations, 
but that different lines respond in different ways, both in amplitude and in 
time scale. These data also showed a complicated relationship between UV 
and X-ray variations and led to the discovery of variable absorption lines in 
the ultraviolet. 

The galaxy NGC 4151 was also monitored spectroscopically in the optical 
at Lick Observatory by Antonucci & CohcnD. They found that the Balmer 
lines seemed to respond to continuum variations on a time scale less than 
around one month (their typical sampling interval). They also found that 
relative to H/3 and Ha, the higher-order Balmer lines and HeiiA4686 varied 
with higher amplitudes. 

Arakelian (Akn) 120 was the first higher- luminosity Seyfert that was mon- 
itored fairly extensively in the opticalE-xE3 as a_res.ult of dramatic Balmer-line 
profile changes that had been detected earlierE3'E3'Ej. Peterson et al.Ej found 
that the time scale for the response of H/3 to continuum variations suggested 
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a BLR size of less than 1 light month across. This was a surprising result as it 
suggested that there was a serious problem with existing estimates of sizes of 
the BLR that were based on photoionization equilibrium modeling, as these 
indicated the BLR should be about an order-of-magnitude larger than this. 
The upper limit on the BLR size was similar to that obtained by Antonucci 
& CohenEI for NGC 4151, but because Akn 120 is a higher-luminosity source, 
the monthly sampled data provided a more critical challenge to BLR models. 

Not surprisingly the results from these earlier monitoring programs were 
controversial. Several observational problems could be identified: 

1. Undersampling of the variations. The variations tended to be undersam- 
pled because the original programs for monitoring Seyfert galaxies were 
designed for BLRs that were thought to be many light months in size. 
For example, in the case of Akn 120, Peterson et al.E3 were looking for 
profile structures that were expected to cross the line profile on a time 
scale of a year or so, and monthly observations should have been sufficient 
to carry out this program. There was certainly at the time no reason to 
believe that higher sampling rates should be required; indeed, proposals 
to observe AGNs as often as once per month were sometimes deemed 
to be "oversampled" by telescope allocation committees! There is no 
obvious algorithm to determine whether or not the variations that have 
been observed are undersampled, but it is quite obvious that if the re- 
sults depend on individual data points, the light curve is almost certainly 
undersampled and any conclusions drawn must be eyed with suspicion. 
A very simple operational criterion for adequate sampling is that if the 
results do not change much when individual points are removed from the 
light curve, the light curve is probably not seriously undersampled. A 
nice simple test is to divide a light curve into two parts, one comprised of 
the even- numbered points (i.e., second, fourth, etc., in the time-ordered 
series) and one comprised of the odd-numbered points. If the two light 
curves are still very similar, i.e., the important features appear in both 
light curves, then the original light curve is probably adequately sampled. 

2. Low S/N of the light curves. If the detected variations are not large 
compared to the signal-to- noise ratio (S/N) of typical data, then spuri- 
ous results can be obtained. Stated another way, F vaI (Eq. (j^)) must 
be significantly greater than zero. This was a serious issue in the case 
of some of the earlier data obtained with image dissector scanners and 
Reticon arrays, for which uncertainties in AGN line and continuum fluxes 
were typically around 8-10%. To a large degree, this problem has been 
obviated by use of CCDs, for which typical errors in the 1-3% range are 
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routinely achieved. 



3. Systematic errors. There are two sometimes-related types of insidious 
errors that can adversely affect time series analysis: (a) correlated con- 
tinuum/line errors, and (b) aperture effects. 

Correlated errors are due to systematic flux-calibration errors. Basically, 
if the flux calibration of a spectrum is incorrect, both the continuum and 
emission-line fluxes measured from it will be in error in the same sense; if 
the calibration is too high, both the emission-line and continuum fluxes 
will be too high. This introduces an artificial correlation between the 
line and continuum at zero lag, and can thus bias the measurement of 
the true lag between them to artificially small values. 

Aperture effects occur when the amount of flux entering a spectrograph 
is not fixed, on account of pointing or guiding errors, or variations in 
seeing in the case of ground-based observations. In point-like sources 
like stars, this affects only the overall photometric accuracy. In nearby 
AGNs, however, both the narrow-line region (NLR) and host-galaxy are 
spatially resolved, and the aperture geometry, centering and guiding, and 
seeing variations can lead to apparent spectral variations. Depending on 
their nature, aperture effects can cause either correlated, uncorrelated, 
or even anticorrelated errors in the line and continuum fluxes. 

2.6 Spectrophotometric Flux Calibration 

In this section, we will outline some of the important considerations for 
flux calibration of spectroscopic monitoring data, directed primarily towards 
ground-based optical observers. 

In AGNs, aperture effects arise because the source is comprised of mul- 
tiple components, some of which have angular structure on scales similar to 
the width of the point-spread function (PSF). The basic requirement for ac- 
curately flux-calibrating AGN spectra is that stable fractions of light from the 
AGN continuum source and BLR (both point-like at even the C.'Ol level) and 
the NLR and the host galaxy (both extended even at arcsecond levels) must 
enter the aperture. We note that even in the case of space-based observa- 
tions, there is a trade-off between aperture size and pointing uncertainty: the 
goal is to minimize the amount of host-galaxy starlight entering the aperture 
(arguing for a smaller aperture), while ensuring that the amount of admitted 
starlight is constant (arguing for an aperture large relative to the pointing ac- 
curacy of the telescope). In the UV, however, this is a less-significant problem 
because the host galaxy is so much fainter than the AGN itself. 
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Mitigation of aperture effects is one of the most important considerations 
for a ground-based monitoring program. This is important to keep in mind, 
as most AGN observers are used to background-limited observations of point 
sources, and this calls for adjusting the slit width for variations in seeing in 
order to optimize the S/N of the data. However, this is exactly the wrong thing 
to do if you are monitoring Seyfert galaxy variations, since calibration accu- 
racy is almost always determined by systematics rather than photon statistics. 
Most observers do know that you need to open up the aperture for absolute 
spectrophotometry, however, and this is precisely what needs to be done in 
monitoring programs, as we will explain below. 

The standard method of absolute flux calibration is to determine from 
observations of spectrophotometric standard stars how counts per second per 
pixel on the detector translates to flux per unit wavelength. However, stan- 
dard absolute spectrophotometry is far too inaccurate for ground-based AGN 
monitoring; the typical accuracy that can be achieved on photometric nights 
is about 5%. Furthermore, at most observing sites, only a relatively small 
fraction of nights are of sufficient quality and stability for absolute spectropho- 
tometry to be useful. Even at a good site like Kitt Peak, fewer than 1/3 of 
the nights are photometric. However, more than 2/3 of all nights (including 
bright time) are good enough to get quality spectra of bright AGNs with a 
1-2-m telescope with a typical photometric accuracy about 15%. Since we 
need a high rate of sampling and do not wish to discard more than half of 
the potentially usable nights, secondary higher-accuracy calibration methods 
must be used. In ground-based AGN monitoring, there are two commonly 
used high-accuracy flux-calibration methods, which we outline below: 

1. Relative spectrophotometry. This method involves simultaneous observa- 
tion of a nearby field star (that must be shown to be non- variable!) in 
long-slit mode. Absolute calibration can be achieved by calibrating spec- 
tra of the comparison star in the usual fashion described above on nights 
of photometric quality. This method has been used to attain internal 
(relative) accuracies of ~ l-2%0J'Ea. The limitations of this method are 
(a) an extremely accurate slit response function (i.e., relative sensitivity 
as a function of position in the slit) is required, (b) a large slit and high 
pointing and guiding accuracy is required to avoid aperture/seeing effects, 
and (c) this method is difficult to use shortward of ~ 4000 A because the 
requirement of fixed position angle (to observe the AGN and comparison 
star simultaneously) leads to atmospheric dispersion problems. 

2. Internal calibration. This is based on using a constant-flux component in 
the spectrum (e.g., narrow emission lines) to refine absolute calibration. 
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this way, internal accuracies of ~ 2% have been routinely achievecH 
. The [Om] AA4959, 5007 doublet is excellent for this method because 
these lines are strong, relatively unblended, and close in wavelength to H/3 
and HenA4686. Unlike the broad lines, the narrow lines are almost al- 
ways expected to be constant in flux. The light-crossing time for the NLR 
is large (typically 100-1000 years) as is the recombination time (about 
100 years; see Eq. ( |l7j ) in Sec. |J), so any short-term variability is smeared 
out. This method also has limitations: (a) most narrow lines are blended 
and/or weak, (b) the method works best at reasonably high spectral res- 
olution, when the narrow-lines are at least partially resolved, and (c) in 
nearby AGNs, the NLR is often extended on scales of arcseconds, which 
introduces the possibility of aperture effects. 

The last of these problems deserves serious consideration. If the NLR is 
extended, measuring a consistent flux is difficult due to aperture effects, as 
described above; either centering and guiding errors or differences in the seeing 
between observations can alter the calibration of the spectrum. However, it 
is possible to understand the nature of these effects so that in principle one 
can correct for them, or at minimum one can pursue mitigation strategies to 
reduce their influence. To do this, one must know (a) the surface brightness 
distributions of the extended components (the host galaxy and the NLR) and 
(b) the PSF, including the seeing component, for the instrument at the time 
of the observations. 

With a good model for the surface brightness distributions, we can still use 
the riarrnw-liryag for flux calibration, and then apply a seeing-dependent flux 
correctioroS. Even if the seeing is not measured, such models can be used 
to give us an idea of how badly our observations might be adversely affected 
by seeing variations. As an example, the corrections to the point-source fluxes 
(broad lines and AGN continuum) and host galaxy flux that need to be applied 
in the case of NGC 4151 are shown in Fig. ||for different aperture geometries. 
We can readily conclude from these diagrams that seeing-dependence of the 
correction factors is smaller for larger apertures, as can be expected intuitively. 
A useful generalization is that the corrections for seeing effects are generally 
not necessary if the aperture in both dimensions is larger than the worst 
seeing, as measured by the full-width at half-maximum (FWHM) of a stellar 
image. 

2. 7 Requirements for a Successful Monitoring Program 

In order to avoid misleading systematic errors of the sort discussed thus far, 
care must be taken in designing an observational monitoring program. The 
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Figure 5. Aperture corrections (as denned by Peterson et al. 68 ) for NGC 4151. These are 
multiplicative adjustments that must be made to (a) the broad emission-line and AGN 
continuum fluxes (top panel) and (b) the host-galaxy starlight contribution (bottom panel) 
as a function of seeing when the spectra are calibrated by assuming a constant flux in 
[Om] AA4959, 5007. The reference point is a 5" X 7'.'5 aperture and seeing FWHM = T.'3. 
For sufficiently large apertures, the correction is negligible under even quite poor observing 
conditions; for small apertures, seeing effects can produce apparent variability. The starlight 
correction factor for the 10" X 13" aperture has a value of ~ 1.6 and varies only weakly 
with seeing. 



fairly obvious basic requirements that we should keep in mind are: 

1. The monitoring program must have sufficient sampling. Specifically, the 
observations should be closely spaced in time relative to physical time 
scales of interest. In terms of total duration, the program should be at 
least three times as long as the longest time scale to be investigated. 

2. The S/N of the light curve must be high enough to detect variations on 
the shortest interesting physical time scale. Inspection of Fig. for ex- 
ample, demonstrates that if you can measure continuum fluxes to only 
10% accuracy, 3er detections of continuum variability in Seyfert 1 galaxies 
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will be extremely rare on time scales as long as months. 

3. Systematic effects must be obviated. The data must be as homogeneous 
as possible. The easiest way to do this is to acquire the data with a single 
instrument in a stable configuration that is suitable for all the observing 
conditions expected over the duration of the program. 



3 Theory of Reverberation Mapping 

After emission^ne-variability was detected, it became clear to a number of 
investigatorsO'c3'l3't2l^ that the kinematics and the geometry of the BLR can 
be tightly constrained by characterizing the emission-line response to contin- 
uum variations. The time delay between continuum and emission-line vari- 
ations are ascribed to light travel-time effects within the BLR; the emission 
lines "echo" or "reverberate" to the continuum changes. The Blandford & 
McKee paperQ, regarded as the seminal paper in the field, first introduced the 
term "reverberation mapping" to describe this prpcess. Reviews of progress-La 
reverberation mapping are provided by Petersonl£3 and Netzer & PetersonEa. 

Because the emission lines vary with the UV/optical continuum (with 
some small time delay), we can immediately draw several important conclu- 
sions: 

1. The line-emitting clouds are close to the continuum source, i.e., the light- 
travel time across the BLR is small. 

2. The line-emitting clouds are optically thick. If they were optically thin, 
the line emission from them would change little as the continuum varied. 

3. The observable UV/optical continuum variations are closely related to 
variations of the ionizing continuum (A < 912 A, hv > 13.6 eV). 

This information provides the crucial underpinning for the assumptions 
detailed below. But before we proceed any further, we need to remind our- 
selves of some of the basic characteristics of the line-emitting gas, much of 
which we derive by comparing the observed line spectra with the predictions 
of photoionization equilibrium computer codes (see the contribution by Netzer 
for a more complete discussion). In general, photoionization equilibrium mod- 
els of the line-emitting clouds are parameterized by the shape of the ionizing 
continuum, elemental abundances, and an "ionization parameter" 

V=^, (13) 
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where 

r°° t 

Q(H) = / -fdu (14) 

is the number of hydrogen-ionizing (hv\ = 13.6 eV) photons emitted per sec- 
ond by the central source. 

Essentially, U characterizes the ionization balance within the cloud, as 
<3(H)/r 2 is proportional to the number of photoionizations occurring per sec- 
ond at the incident face of the cloud, and n c is proportional to the recombi- 
nation rate. For a typical Seyfert 1 galaxy (NGC 5548) 

Q(H) re 10 54 h^ 2 photons s" 1 , (15) 

where ho is the Hubble constant in units of 100 km s _1 Mpc -1 . By taking the 
electron density to be n c ps 10 11 cm _3 and the distance to the central source 
to be 10 light days, we find U ~ 0.1. The line-emitting clouds are ionized to 
the "Stromgren depth" , 

TTr 

R s = » 10 11 cm , (16) 

aB"c 

at which point virtually all of the incident ionizing photons have been ab- 
sorbed. In BLR clouds, the region beyond the Stromgren zone is partially 
ionized, n(H + ) ps 0.1n(H°), and heated primarily by far infrared radiatiorO 
and partly by high-energy X-raysc3. In Eq. (|l6|), as is the Menzel-Baker case 
B recombination coefficicntEa 

3.1 Reverberation Mapping Assumptions 

There are a number of simplifying assumptions that we can now make: 

1. The continuum originates in a single central source. The size of an accre- 
tion disk around a supermassive (say, 10 7-8 M Q ) black hole is of the order 
of 10 13_14 cm. A typical size for the BLR in the same system would be 
of order a few light days, i.e., ~ 10 16 cm. Note that it is specifically not 
required that the continuum source emits radiation isotropically, though 
this is a useful starting point. The "point-source" assumption greatly 
simplifies the reverberation process. However, we should mention that 
the point-source assumption is probably not applicable to X-ray rever- 
beration, as the Fe Ka emission and the hard X-rays that drive this line 
probably arise in regions that are of similar size, and perhaps co-spatiaO. 

2. Light-travel time is the most important time scale. We assume specifi- 
cally that emission-line clouds respond instantaneously to changes in the 
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continuum flux. The time scale to re-establish photoionization equilib- 
rium is the recombination time, 

W = Kaer 1 « 40 ( ^^-3 ) s. (17) 

The time it takes a Lyman a photon to diffuse outward Jrom the 
Stromgren depth is about 20 times the direct light-travel timc£3, 

r diff = 20^ « 20 — « 60 ( - " c 3 ) ~* s. (18) 
c n c a B VlO^cm - - 3 / 

We also need to carry out our reverberation-mapping experiment on a 
time scale short enough that the structure of the BLR can be assumed 
to be stable. The dynamical, or cloud-crossing, time scale for the BLR 
is typically 

T 

Tdyn = 77 ~ 3 - 5 years, (19) 

'FWHM 

where Vf WHM is the Doppler width of the broad line for which the response 
time t = r/c has been measured. Any reverberation-mapping experiment 
has to be short relative to the crossing time or the structural information 
might be washed out by cloud motions. 

3. There is a simple, though not necessarily linear, relationship between the 
observed continuum and the ionizing continuum. 

3.2 The Transfer Equation 

Under these assumptions, the relationship between the continuum and emis- 
sion lines can be written in terms of the "transfer equation" , 

/oo 
*(V z> r)C(t-T)dT . (20) 
-00 

Here C(t) is the continuum light curve, L(V z ,t) is the emission-line flux at 
line-of-sight velocity V z and time t, and ^(V z ,t) is the "transfer function" 
at V z and time delay r. Inspection of the transfer equation shows that the 
transfer function is simply the time-smeared emission-line response to a 5- 
function outburst in the continuum. Note that causality requires that the 
lower limit on the integral is r = 0. 

Solution of the transfer equation to obtain the transfer function is a classi- 
cal inversion problem in theoretical physics: the transfer function is essentially 
the Green's function for the system. Unfortunately, it is difficult to find a sta- 
ble solution to such equations, especially when the data are noisy and sparse, 
as they usually are in astronomical applications. 
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In practice, most analyses have concentrated on solving the velocity-in- 
dependent (or 1-d) transfer equation, 



where both ^(r) and L(t) represent integrals over the emission-line width. 

The transfer equation is a linear equation. In reality, however, the re- 
lationship between the observed continuum and the emission-line response 
is likely to be nonlinear. We therefore approximate C{t) = C + AC and 
L{t) — L + AL, where C and L represent constants, usually the mean value of 
the continuum and line flux, respectively. We can then treat deviations from 
the mean as linear perturbations and the equation we actually solve is 



Most existing data sets are inadequate for transfer-function solution, and 
simpler analyses are used. The most commonly used tool in analysis of AGN 
variability is cross-correlation, which will be discussed in detail in Sec. ||. 

3.3 Isodelay Surfaces 

Suppose for the moment that the BLR consists of clouds in a thin spherical 
shell of radius r. Further suppose that the continuum light curve is a simple 
8- function outburst. Continuum photons stream radially outward and after 
travel time r/c, about 10% of these photons (using a typical "covering factor") 
are intercepted by BLR clouds and are reprocessed into emission-line photons. 
An observer at the central source will see the emission-line response from 
the entire BLR at a single instant with a time delay of 2r/c following the 
continuum outburst. At any other location, however, the summed light-travel 
time from central source to line-emitting cloud to observer will be different 
for each part of the BLR. In the case of a 5- function outburst, at any given 
instant, the parts of the BLR that the observer will see responding are all 
those for which this total path length is identical; at any given time delay, 
the part of the BLR that the observer sees responding is the intersection of 
the BLR distribution and an "isodelay surface." Astronomers, on account of 
their familiarity with conic sections, can readily recognize that the shape of 
the isodelay surface is an ellipsoid with the continuum source at one focus and 
the observer at the other; the light-travel time from central source to BLR 
cloud to observer is constant for all points on the ellipsoid. Since the observer 
is virtually infinitely distant from the source, the isodelay surface becomes a 
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(22) 




Figure 6. The circle represents a cross-section of a shell containing emission-line clouds. 
The continuum source is a point at the center of the shell. Following a continuum outburst, 
at any given time the observer far to the left sees the response of clouds along a surface 
of constant time delay, or isodelay surface. Here we show five isodelay surfaces, each one 
labeled with the time delay (in units of the shell radius r) we would observe relative to the 
continuum source. Points along the line of sight to the observer are seen to respond with 
zero time delay. The farthest point on the shell responds with a time delay 2r/c. 



paraboloid, as shown schematically in Fig. |6| The figure shows the BLR as 
a ring intersected by several isodelay surfaces, labeled in terms of their time 
delay in units of r/c. Relative to the continuum, points along the line of sight 
to the observer are not time delayed (i.e., t = 0). Points on the far side of the 
BLR are delayed by as much as 2r/c, the time it takes continuum photons 
to reach the BLR plus the time it takes line photons emitted towards the 
observer to return to the central source on their way to the observer. 

Essentially, the transfer function measures the amount of line emission 
emitted at a given Dopplcr shift in the direction of the observer as a function 
of time delay r. The value of the transfer function at time delay r is computed 
by summing the emission in the direction of the observer at the intersection 
of the BLR and the appropriate isodelay surface. For a thin spherical shell, 
the intersection of the BLR and an isodelay surface is a ring of radius r sin#, 
where the polar angle 9 is measured from the observer's line of sight to the 
central source, as shown in Fig. |?]. The time delay for a particular isodelay 
surface is the equation for an ellipse in polar coordinates, 

t = (1 + cos 9)r/c , (23) 

as is obvious from inspection of Fig. [?|. The surface area of the ring of radius 
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Figure 7. The upper diagram shows a ring (or cross-section of a thin shell) that contains 
line-emitting clouds, as in Fig. 6l An isodelay surface for an arbitrary time is given; the 
intersection of this surface and the ring shows the clouds that are observed to be responding 
at this particular time. The dotted line shows the additional light-travel time, relative to 
light from the continuurnspurce, that signals reprocessed by the cloud into emission-line 
photons will incur (Eq. (E3)). In the lower diagram, we project the ring of clouds onto 
the line-of-sight velocity/time-delay (Vz,t) plane, assuming that the emission-line clouds 
in the upper diagram are orbiting in a clockwise direction (so that the cloud represented 
by a filled circle is blueshifted and the cloud represented by an open circle is redshifted). 



rsin# and angular width r dO is 2ir(r sin#) r d9, and assuming that the line 
response per unit area on the spherical BLR has a constant value Eq, the 
response of the ring can be written as 

= 2we r 2 sin 6 d9 , (24) 

where < 9 < 2ir. From Eq. (p3|), we can write 

dr = -(r/c) sin 9 d9 , (25) 
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so putting the response in terms of r rather than 9, we obtain 



*(r)dr = m(9) 



(10 



dr = 2irerc dr (26) 



for values from r = (0 = 27r) to r = 2r/c (0 = 0). The transfer function for 
a thin spherical shell is thus constant over the range < r < 2r/c. 

3-4 Transfer Functions for a Variety of Simple Models 

The simple analytic calculation in the last section was intended to be illustra- 
tive and serve as a reference point. We will now expand on this with more gen- 
eral geometries, and also incorporate information from Doppler motion along 
the line of sight. We will derive some transfer functions for other simple mod- 
els, focusing on three: (a) systems of clouds in circular Keplerian orbits, illumi- 
nated by an isotropic continuum, (b) biconical outflows, and (c) disks of ran- 
dom inclination. All these are physically plausible, and can produce "double- 
peaked" emission-line profiles, which are sometimes seen in AGNs, though not 
all of these models necessarily do this. We will start with the simplest models 
and progress to more complicated models. Much of this discussion is drawn 
from discussions of transfer functions in the literatureEj'EEl'Ej'LJ'EJ. 

Suppose line-emitting clouds are on a circular orbit at inclination i — 90°; 
imagine that the circle in Fig. 0a represents this orbit seen face on. The line 
response from the clouds at the intersection of an arbitrary isodelay surface 
and the circular orbit will be at time delay r = (1 + cos 9)r/c and line-of-sight 
velocities V z = ±Vo r bSin#, where Vorb = (GM/r) 1 / 2 , the circular orbital 
speed. It is easy to see that the circular orbit projects to an ellipse in the line- 
of-sight velocity/time-delay (V z , r) plane with semiaxes Vorb and r/c, as shown 
in Fig. ^b. This simple example is important because it is straightforward to 
generalize it to both disks (rings of different radii) and shells (rings at different 
inclinations). 

First we consider the generalization to a shell. We can construct a shell 
from a distribution of circular orbits, with inclinations ranging from i = 0° 
to i = 180°. As we decrease the inclination of the circular orbit in Fig. ^a 
from i — 90°, we see that the range of time delays will decrease from [0, 
2r/c] to [(1 — sini)r/c, (1 + srni)r/c], and similarly the line-of-sight velocity 
range will decrease from [— V OI b, +V OI b] to [—Vorb sin j, -I- Vorb sin?], as we show 
schematically in Fig. |8[ At the limiting case i = 0°, the time delays all contract 
to r/c, since the light travel-time paths for all points on a face-on ring are the 
same, and the velocities all contract to V z = 0, because the orbital velocities 
are now perpendicular to the line of sight. Thus, the transfer function looks 
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Figure 8. Here we take the circular ring shown in Fig. [7| and show how its projection on 
V z ,t changes as the inclination is decreased from 90°; in the time-delay direction, the axis 
contracts to 2rsini/c, and in the velocity direction it decreases to 2V or b sini. 



like a series of ellipses as in Fig. ||b that with decreasing inclination contract 
down to a single point at V z = and t = r/c when i = 0°. We can construct 
such a transfer function by using a Monte Carlo method that places BLR 
clouds randomly across the surface of the shell, and the result we get is shown 
in Fig. |^, which corresponds to a thin shell of radius 10 light days and a central 
mass of 10 8 M©. We have also integrated this transfer function over r and V z 
to obtain the emission-line profile and the one-dimensional transfer function, 
respectively. In the particular case of a thin spherical shell, we see that both 
of these are simple rectangular functions, as we showed analytically for the 
one-dimensional transfer function (Eq. (p6|)). 

The first complication that we might consider is anisotropic line emission 
by the BLR clouds. Physically, this will occur if the BLR clouds themselves 
are optically thick in the lines as well as the continuum. In this case, most 
of the line radiation emitted by the clouds will be from the side of the cloud 
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Figure 9. Transfer function for a thin spherical shell. The upper left panel shows in grey 
scale the two-dimensional transfer function, i.e., the observed emission-line response as a 
function of line-of-sight velocity V z and time delay r. The upper right panel shows the one- 
dimensional transfer function, i.e., the two-dimensional transfer function integrated over 
V z , which is the response of the total emission line as a function of time (cf. Eq. (p6[)). 
The lower left panel shows the emission-line response integrated over time delay; this is the 
profile of the variable part of the line. The specific model here is comprised of emission-line 
clouds in circular orbits of radius r = 10 light days and random inclinations (0° < i < 180°), 
orbiting around a central black hole of mass 10 s Mq. In this model, the line emission from 
each cloud is isotropic (A = 0). 



facing the continuum source, i.e., the line emission is directed back towards 
the continuum source. A simple parameterization of asymmetric line emission 
is to describe the apparent emissivity of a cloud as 

e(6) = £ (1 + Acos9) , (27) 

where Eq is constant and the parameter ^4 = for isotropic emission and A = 1 
for completely anisotropic emission; the latter case is appropriate for spherical 
clouds with inward- facing surfaces that are uniformly bright. In principle, A 
can be estimated by photoionization modeling, though in practice the values 
are highly uncertain on account of limitations in the accuracy of the radiative 
transfer codes (see the contribution by Netzer). It is certainly expected that 
Ami for Lya, at least, and models suggest approximate values 4 k 0.7 
for Civ and A ss 0.8 for H/3o. The main effect of anisotropic emission is to 
increase the measured lag for a given geometry because the apparent response 
of the near side of the BLR is suppressed. For a thin shell of radius r, the 
mean time delay is r = (1 + A/3)r/c. Fig. [To] shows the transfer function 
for the same thin spherical shell model of FigrTI, but in this case with highly 
anisotropic (^4 = 1) line response. 

In addition to anisotropic line response, we can also consider anisotropic 
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Figure 10. Transfer function for a thin spherical shell, as in Fig. f), except with completely 
anisotropic (A = 1) line emission, i.e., all of the emission line flux from each cloud is directed 
back towards the continuum source. 



illumination of the BLR clouds by the continuum source. As an example, 
consider the case where BLR clouds are illuminated by biconical beams of 
semi-opening angle to; the limiting case as w approaches zero would be a 
narrow jet-like pencil beam, and the case u> = 90° corresponds to an isotropic 
continuum. We start by examining the response of an edge-on (i = 0°) circular 
ring, as we show in Fig. [H], which is exactly like Fig. |7[ but with only part 
of each orbit illuminated by the continuum source. We earlier generalized the 
result for a shell, going from Fig. [?] to Fig. || as we decreased the inclination 
of the ring; doing this again, we see in Fig. ^ how the two-dimensional 
transfer function is altered by anisotropic illumination. We note specifically 
the absence of any response near r = and t = 2r/c since the bicones shown 
do not illuminate BLR clouds directly along the observer's axis. Similarly 
there is also no response around r = r/c due to the absence of response of 
clouds at w 90°. 

So far we have restricted our attention to "thin" geometries, i.e., single 
orbits and shells. Generalization to "thick" geometries, e.g., disks and shells 
with different inner and outer radii, is trivial: the response of a disk can be 
computed by integrating over a series of circular orbits, and the response of 
a thick shell is obtained by integrating over a series of shells. In Fig. we 
illustrate this concept by showing the response from two rings. Thus far, the 
free parameters we have dealt with are the radius r, the line anisotropy factor 
A, and in the case of non-spherically symmetric systems, the inclination i and, 
in the case of biconical illumination, the beam half-angle lo. Thick geometries 
now introduce inner and outer radii, r m and r out , respectively, and a distance- 
dependent responsivity per unit volume (or per unit area, for a disk), which 
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Figure 11. A ring-like line-emitting region and its projection into the V z ,t plane is shown, 
as in Figs. [?j and B though in this case with a biconical continuum geometry. In the case 
shown here, the beam semi-opening angle is uj = 30°, and the beam is inclined to the 
observer by i = 45°. Again assuming that clouds on the ring are orbiting in the clockwise 
direction, we highlight in the lower diagram the regions in the (V z , r) plane that will produce 
an emission-line response; there is no emission-line response from clouds not illuminated by 
the continuum beam. 



we can parameterize as e(r) oc eqt p . The index (3 allows us to condense 
model-dependence into a single parameter that will account for effects due 
to geometrical r~ 2 dilution of the continuum, a distance-dependent covering 
factor, etc. In Figs. [l4| and [if] we show transfer functions for thick spherical 
shell systems with A = and identical values of r- m and r out , but differing 
radial-response indices /?; the effect of increasing (3 is to enhance the relative 
response of material at larger radii; the limiting cases where (3 — > — oo and 
(3 — ► +oo correspond to the response functions of thin shells of radius r m and 
''out, respectively. We will show additional examples below. 

We consider now a thick shell system that is illuminated by an anisotropic 
continuum source. Again, we assume that the line-emitting clouds are in 
circular Keplerian orbits of random inclination. We show examples that are 
identical except for the continuum beam width and inclination in Figs. [l6| and 
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Figure 12. Transfer function for a thin spherical shell, as in Fig. g, except with an anisotropic 
continuum source. In the example shown, the beam opening angle is ui = 30° and the 
inclination is i = 45°, as in Fig. |llj. 




Line-of-sightvelocity V 



Figure 13. Two ring-like regions (as in Fig. and their projected response on the (V z ,r) 
plane. 



p7| . An important thing to notice is that in one case (Fig. |T^) the observer is 
outside of the continuum beam (i.e., i > uj) and in the other case (Fig. |T7| ) 
the observer is inside of the continuum beam (i.e., i < uj); when the observer 
is inside the beam, the line profile is single-peaked, and when the observer is 
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Fi gure 14. Transfer function for a thick spherical shell with v in — 2 It-days, r ou t = 10 
lt-days, and radial responsivity index j3 = 0. 



Figure 15. Transfer function for a thick spherical shell exactly as described in Fig. |^|, 
except with radial responsivity index j3 = —4. 



outside the beam, the line profile is double-peaked. 

Another simple kinematic model for the BLR consists of clouds in spher- 
ical outflow. The transfer functions for such cases are quite distinctive; an 
example of a two-dimensional tranfer function for a kinematic field with ra- 
dial velocity V r oc r for r less than some maximum distance r out is shown 



in Fig. 18. This velocity field corresponds to either a ballistic outflow or a 
flow undergoing constant acceleration. At t = we see response from all the 
material along the line of sight to the continuum source, which ranges from 
V z = for the material closest to the central source to V z = — V(r out ) for 
the gas farthest from the continuum source. As r increases, we begin to see 
response from the far side of the BLR where the line of sight velocities V z 
become positive. The range of line-of-sight velocities decreases as the isode- 
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Figure 16. Transfer function for a thick spherical shell with rj n = 2 lt-days, r ou t = 10 
lt-days, radial responsivity index j3 = —2, and isotropic line response A = 0. In this model, 
the shell is illuminated by a biconical continuum with semi-opening angle u) = 30° and 
inclination i = 45°, i.e., the observer is outside the continuum beam. Note the double- 
peaked line profile, and contrast this with the model shown in Fig. llTj. 



Figure 17. Transfer function for a thick spherical shell exactly as described in Fig. [If], except 
that in this model the shell is illuminated by a biconical continuum with semi-opening angle 
u> = 75° and inclination i = 15°, i.e., the observer is well inside the continuum beam. Note 
the single-peaked line profile, and contrast this with the model shown in Fig. hd 



lay surfaces get farther from the line of sight. At t — r out /2, the isodelay 
surface no longer crosses any clouds moving towards the observer, and the gas 
moving fastest away from the observer is that along the line of sight (8 = 0) 
with V z — V(r ut)/2. At r = 2r out /c, only the response from the antipodal 
point is seen, and the transfer function is contracted to the single point at 
[2r out /c,V(r out )}. 

The case of biconical outflows (which are possibly relevant, as they are 
certainly seen in the NLR and might well apply to at least a component of 
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the BLR) can be dealt with by restricting the response to certain values of 9; 
an example is shown in Fig. |l9| . 

We have now seen that both orbital and outflow models can produce 
similar emission-line profiles; if the line-emitting gas is confined to a bicone, 
either because of the gas distribution or the ionizing-photon distribution, one 
can get a single-peaked or double-peaked line profile. The two situations can 
be easily distinguished, however, by their two-dimensional transfer functions 
(or equivalently, the combination of their one-dimensional transfer functions 
and their line profiles). In Fig. [2^, we directly compare the one-dimensional 
transfer functions and line profiles for two thick-shell models: (1) emission- 
line clouds in a biconical outflow and (2) clouds in circular Keplerian orbits 
of random inclination, illuminated by a biconical continuum source. In both 
cases, the beams (one radiation, one matter) have the same half-opening angle 
(uj = 40°) and two different inclinations are shown; i = 25° in the top row (i.e. 
the observer is in the beam, as indicated in the left-hand column illustrations 
of the geometry), and i — 65° in the bottom row (i.e., the observer is out of 
the beam). The distribution of line-emitting clouds is the same, regardless of 
how the clouds are moving, so in these two cases the one-dimensional transfer 
functions ought to be very similar, which is indeed the case, as seen in the 
middle column of Fig. |2^. The right-hand column shows the line profiles, 
which are however very different. Consider the case i = 25° (top row): in 
the case of outflow, the line-emitting material in the beam is moving radi- 
ally outward, giving relatively highly blueshifted (near side) and redshifted 
(far side) emission, but little emission near V z = 0, since there is no line- 
emitting material moving transverse to the line of sight. In the case of clouds 
in circular orbits illuminated by an anisotropic beam, the cloud motions are 
perpendicular to their radial vectors, so most of the line-emitting material is 
at low Doppler shift as the gas motions through the beams are predominantly 
transverse. Now consider the higher-inclination case (i = 65°; bottom row): 
in the case of radial outflow, the gas motions in this case are now primarily 
transverse so the Doppler shifts are smaller. However, in the case of circular 
orbits, the material in the beam is moving primarily along the line of sight, 
and there is a deficiency of material at V z « 0. Note that either kinematic 
field can give either double-peaked or single-peaked profiles: a simple gener- 
alization is that double-peaked profiles are produced in outflow models when 
the observer's line-of-sight is in the beam (i.e., i^w) and in orbital models 
when the line of sight is out of the beam (i £ u>). Neither the profiles nor the 
one-dimensional transfer functions individually tell us much about the BLR 
kinematics and velocity field, but together they can tell us a lot. Information 
on both r and V z , i.e., the two-dimensional transfer function, greatly reduces 



in The Starburst-AGN Connection 2001, (Singapore: World Scientific) 32 



Figure 18. Transfer function for a spherical outflow, with outflow velocity V oc r. For this 
model, r in = 0.1 lt-days, r out = 10 lt-days, V(r out ) = 10000km s _1 . 



Figure 19. Transfer function for a biconical outflow, with parameters as in Fig. [is] except 
that the outflowing gas is confined to a bicone of half- width u> = 30° at inclination i = 45° . 



the ambiguities. 

Finally, for completeness, we consider the case of an inclined disk, as this is 
the classic geometry for producing a double-peaked line profile. The response 
of a disk-like BLR can be computed by integrating the response of rings of 
different radii. In Fig. we show the transfer function and line profile for 
a disk at intermediate inclination (i = 45°). Identical line profiles can be 
obtained for different inclinations simply by suitably adjusting the central 
mass, but the transfer function allows us to reduce the ambiguity between 
possible disk models. 
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Figure 20. Response of two biconical models with different kinematics. The left hand 
column shows the two geometries considered, both with beam half-opening angles ui = 40°; 
the top row is at inclination i = 25° (i.e., the observer is in the beam) and the bottom row is 
i = 65° (i.e., the observer is out of the beam). One-dimensional transfer functions (middle 
column) and line profiles (right column) are shown. The solid line is a biconical outflow 
model, and the dashed line represents circular Kcplcrian orbits of random inclination. The 
one-dimensional transfer functions are identical, but the profiles are greatly different, as 
described in the text. 




Figure 21. Transfer function for a thin disk at intermediate inclination (i = 45°). Note 
that the line profile is double peaked, and the one-dimensional transfer function is similar 
to those for spherical shells. 
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3.5 Transfer Function Recovery 



With data of sufficient quality, transfer functions for different plausible models 
are distinguishable from one another. This is of course, why we place so much 
emphasis on them: if we can determine the transfer function for a particular 
line, we have very strong constraints on the geometry and kinematics of the 
responding region, and if the BLR has a spherically or azimuthally symmetric 
structure, it may be possible to determine the BLR kinematics and structure 
with fairly high accuracy and strongly constrain the emission-line reprocessing 
physics. The operational goal of reverberation mapping experiments is thus 
to determine the transfer functions for various emission lines in AGN spectra. 

How can we determine the transfer functions from observational data? 
Inspection of the transfer equation (Eq. (^0|)) immediately suggests Fourier 
inversion (i.e., the method of Fourier quotients), which was the formulation 
outlined by Blandford & McKeeQ. We define the Fourier transform of the 
continuum light curve C(t) as 

/+oo 
C(t) e-^ft dt (28) 
-oo 

and similarly define n the Fourier transform of the line light curve. By the 
convolution theorems, the transfer function becomes 

L* = x C* , (29) 

so that \P* = L* /C*, and the transfer function is obtained from 

/+oo 
q*(f)e t2 * ft dt . (30) 
-oo 

In practice, however, Fourier methods are not used as they are inadequate 
when applied to data that are relatively noisy (i.e., flux uncertainties are only 
a factor of a few to several smaller than the variations) and which are limited 
in terms of both sampling rate, which is in any case usually irregular, and 
duration. Simpler methods, like cross-correlation (Sec. ^), can be applied 
with success, though with limited results. Cross-correlation, we will sec, can 
give the first moment of the transfer function, but little else. 

In principle, more powerful methods can be used to attempt to recover 
transfer functions.. The most commonly used is the Maximum Entropy 
Method (MEM)Ea. MEM is a generalized version of maximum likelihood 
methods, such as least-squares. The difference is that in the method of least 
squares, a parameterized model is fitted to the data, whereas MEM finds the 
"simplest" (maximum entropy) solution, balancing model simplicity and real- 
ism. Examples of MEM solutions will be shown in Sec. ||. Other methodologies 
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Figure 22. The left-hand panel shows simultaneous measurements of H/3 emission-line and 
optical continuum fluxes for the Seyfert galaxy Mrk 335 69 . These consist of 24 measure- 
ments made on an approximately weekly basis over one observing season. The right-hand 
panel shows the same emission-line fluxes paired with continuum values from 15.6 days 
earlier; i.e., the emission-line fluxes are better correlated with earlier rather than current 
continuum values on account of light travel-time delays. When the time lag is taken into 
account, the fit improves. 



that have been employed for transfer-function solution include suhtractively 
optimized local averages (SOLA)LD and regularized linear inversions. 

4 Cross-Correlation Analysis 

4-1 Fundamentals 

Cross-correlation analysis is the tool most commonly used in the analysis 
of multiple time series. Because its application to astronomical time se- 
ries is often misunderstood and has historically been rather contentious, 
it merits special attention. Important steps in the development of cross- 
correlatiogi^al^g ^aj applied to AGN variability studies can be found in the 

Cross-correlation analysis is basically a generalization of standard linear 
correlation analysis, which provides us with a good place to start. Suppose we 
obtain repeated spectra of one of the brighter Seyfert galaxies, and we want 
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to determine whether or not the variations in the H/3 emission line and the 
optical continuum are correlated (which was an interesting question 20 years 
ago, even before emission-line time delays were considered). The first thing 
you would do is plot the H/3 flux against the continuum flux, as in the left-hand 
panel of Fig. ^2|, which shows that the two variables are indeed correlated. A 
measure of the strength of the correlation is given by the correlation coefficent, 

„_ E^i fa - g) fa - v) ( 31 j 

Eili {?% - xf^j (Jy£=i (vi - y) 2 ) 

where there are N pairs of values (xj, Hi) and their respective means are x and 
y. When the two variables x and y are perfectly correlated, r = 1. If they 
are perfectly anticorrelated, r — — 1. If they are completely uncorrelated, 
r = 0. For the data shown in the left panel of Fig. [2^, r = 0.596; for 24 
pairs of points, as shown here, this means that the correlation is significant 
at the £ 99.8% confidence level (i.e., the chance that the two variables are 
in fact completely uncorrelated and the correlation we find is spurious is less 
than 0.02%. Confidence levels for linear correlation can be found in standard 
statistical tablesQ) 

While this is quite a good correlation, we see something more remarkable 
if we plot both variables as functions of time (i.e., as light curves), as seen 



in Fig. 23, We see that the patterns of variation are very similar, except 
that the emission-line light curve is delayed in time, or "lagged," relative 
to the continuum light curve. It is obvious that the correlation between the 
continuum and emission- line fluxes would be even better if we allowed a linear 
shift in time between the two light curves in order to line up their prominent 
maxima and minima. This is what cross-correlation does. 

The first operational problem in computing a cross-correlation is also 
immediately apparent: since each point in one light curve must be paired with 
a point in the other light curve, it is obvious that the data should be regularly 
spaced. The cross-correlation is then evaluated as a function of the spacing 
between the interval between data points At using the pairs [x(ti), yfc+NAt)] 
for all integers N. Unfortunately, regularly sampled data are almost never 
found in Astronomy; ground-based programs have weather to contend with, 
and even satellite-based observations are almost never regularly spaced in 
time. The essence of the cross-correlation problem in Astronomy is dealing 
with time series that are not evenly sampled. Moreover, the light curves are 
often limited in extent and are noisy. 



For well-sampled series as in Fig. 23, the sampling problem can be dealt 



with in a straightforward fashion. The simple, effective solution is to interpo- 
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Figurp 23. The H/3 emission-line and optical continuum fluxes for Mrk 335, as shown in 
Fig. ba, are plotted as a function of time. It is clear from the figure that the continuum 
and emission-line fluxes are well-correlated, and that the correlation can be improved by a 
linear shift in time of one time series relative to the other. The optimum linear correlation 
occurs by shifting the emission-line light curve backwards by 15.6 days. 



late one series between the actual data points, and use the interpolated points 
in the cross-correlation. We illustrate this schematically in Fig. We can 
then compute the cross-correlation function CCF(r), as shown in Fig. and 
the step size we use for r is now somewhat arbitrary. At each value of the 
lag r we compute r as in Eq. (|3l|). For the example we have been using, we 
find that the CCF is maximized when points in the continuum light curve are 
matched to those in the emission-line light curve with a delay of 15.6 days. 
If we plot the shifted emission-line values versus the continuum values (as 
we have done in the right-hand panel of Fig. p2| ) and again perform a linear 
correlation analysis, we find that the fit has improved, with r — 0.849 and 
xl = 1-76. 



4-2 Relationship to the Transfer Function 

It is useful at this point to examine how the CCF is related to more fun- 
damental quantities. Defined as a function of the continuous variable r, the 
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Figure 24. Continuum and emission-line light curves for Mrk 335, as in Fig. |23|. This 
illustrates the interpolation method commonly used in cross-correlation. In this figure, 
the emission-line light curve is made continuous through linear interpolation between data 
points. Actual continuum observations are then paired with interpolated emission-line 
values to compute the correlation coefficient for a particular time delay. In this example, 
we show interpolated emission-line fluxes that are time-delayed relative to the continuum by 
15.6 days, which is the value at which the cross-correlation function peaks. As a visual aid, 
dotted lines join a few of the data pairs. Notice how the first few points of the emission-line 
series and the last few points of the continuum series remain unused. 



formal definition of the cross-correlation function is 

/oo 
L(t)C(t -r)dt . (32) 
-OC 

If we use Eq. ( ^T| ) to replace L(t) in this equation, we obtain 
CCF(t) = / C(t - t) / *(r')C(t - t 1 ) dr' dt 



/OO pOO 
*(r ) / C{t- T ')C{t-T) dt dr' . (33) 
- oo J —oo 

Comparing the inner integral with Eq. (|32"|), we see that it is the cross- 
correlation of the continuum light curve with itself, i.e., the "autocorrelation 
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Figure 25. The interpolation cross-correlation function for the Mrk 335 data shown in the 
previous figures. The shaded area indicates points with values r > O.SYpeak, which are the 
points used in computing the centroid, which is also indicated. 



function" , 

/oo 
C(t)C(t - r) dt . (34) 
-OO 

Thus, we see that the cross-correlation function is the convolution of the 
transfer function with the continuum autocorrelation function, 

/oo 
*(t')ACF(t - T ') dr' . (35) 
-oo 

The centroid of the CCF is a usually quoted quantity, 

_ Jt CCF(t) dr 

Tcent- J CCF(r) dT ' ( 36 ) 

which is related to the centroid of the transfer function 

JVtt(7-)dr 

T -nt" ^{r)dr ■ [6() 

It is important to note that these two quantities are the same only in the 
limit where the light curves extend to infinity. Operationally, the centroid is 
computed using only points around the most significant peak, usually those 
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points for which CCF(r) > 0.8r pca k, where r pea k is the maximum value of the 
CCF. Sometimes the location of the peak r pC ak (CCF(r pGa k) = r pca k) is the 
statistical quantity used. 

Cross-correlation does not necessarily give a clean and unambiguous mea- 
sure of the relationship between two time series. In particular, the CCF is 
a convolution of the transfer function with the continuum ACF, so the CCF 
depends on the continuum behavior. To illustrate this, we show in Fig. ^6] 
examples of light curves for various power-law PDSs (Eq. (j^)) and the very 
different ACFs computed for each example. 



Figure 26. The left-hand column shows artificial light curves generated with different power- 
law indices a for the power-density spectra, P(f) oc f~ a . All light curves consist of 2048 
points, normalized by setting -F V ar = 0.3 (Eq. (Q)). The right-hand column shows the 
autocorrelation functions for these light curves. Cases shown, top to bottom, are: (a) 
a = 0. There is no correlation from one point to the next, as seen clearly in the ACF, 
which is essentially zero except at r = 0. This is known as "white noise." (b) a = 1. 
This is the power spectrum characteristic of a random walk; each continuum point now has 
some dependence on history. This case is sometimes called "flicker noise" (see Press 78 for 
further interesting discussion of this case), (c) a = 1.5. This appears to be rather typical 
of AGN variability, (d) a = 2. This power spectrum is characteristic of a supcrpositioning 
of independent events, which is often referred to as "shot noise." 



4-3 Discrete Correlation Methods 

There are some circumstances under which one might not be able to reason- 
ably interpolate between gaps in data. This can occur (a) when there are a 
few large gaps in otherwise well-sampled series or (b) when there is reason 
to believe that the variations might be at least somewhat undersampled. In 
these cases, interpolation might be highly misleading, and another methods 
ology needs to be employed. The "discrete correlation function" (DCF)EJ 
method is one where no assumption about light curve behavior needs to be 
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Figure 27. Part of the light curves from the previous figures, expanded to show DCF bins. 
The particular example shown is for a time-lag bin width At = 4 days and shows the 
location of the bins, pairing real continuum points with emission-line points within a 4-day 
window, centered at a shift of 4 days (i.e., emission-line data points shifted by 2-6 days from 
the continuum points). Note that in this expanded region there are only three emission- 
line points that fall into the bins; in the range shown here, only two continuum and three 
emission-line points contribute to the computation of the correlation coefficient. 



made. The DCF method deals with irregularly sampled data by binning the 
data in time, as illustrated schematically in Fig. This is an alternative 
approach to the irregular sampling requirement: instead of requiring that 
points contributing to CCF(t) are separated in time by exactly the interval 
r, we time-bin the data by pairing points with time separations in the range 
r ± Ar/2, where At is the width of one time bin. Choice of the binning 
window is a free parameter, and two examples are shown in Fig. |28|. 

The principal virtues of the DCF method are (a) that only actual data 
points are used and (b) that it is possible to assign a statistical uncertainty 
to the value of the correlation coefficient in each bin. The relative weakness 
of the DCF is that the data are in some ways underutilized, as is evident in 
Fig. ^7]; for a small data set, the DCF method might completely miss a real 
correlation, although-ii is less like to find a spurious correlation than is the 
interpolation methooO. 

One difficulty of the DCF method is that the number of points per time 
bin can vary greatly, as can be easily inferred from inspection of Figs. ^7] and 
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Figure 28. Cross-correlation functions for the Mrk 335 light curve shown in Fig. 23l The 
DCF values are shown as points with error bars, and the interpolation CCF, as in Fig. Gq, 
is shown as a solid line. The upper panel shows the DCF with a bin width of At = 4 days, 
and the lower panel shows a bin width of At = 8 days; with At = 4 days, some bins have 
no data (those with values set to r = —1), but with Ar = 8 days, the DCF is somewhat 
underresolved. 



. One solution to this is to vary the width of the time bins to ensure that 
there are a statistically meaningful number of points in each bin. A method 
for accomplishing this is the "Z-transformed DCF (ZDCF)"ll 

4-4 Computational Issues 

In the foregoing discussion, we have concentrated primarily on conceptual 
issues without paying much attention to a number pfpCippiputational issues 
that must be dealt with in actually computing CCFsEaej'Ell. 
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1. Edge effects. Only at zero lag do all the points of the time series enter 
into the calculation of the correlation coefficient; at any other lag, points 
at the ends of the time series drop out of the calculation since there are no 
points in the other series to which they can be matched (see Fig. This 
means that at larger and larger lags, fewer and fewer points contribute 
to the calculation of r. This has two significant consequences: 

(a) Normalization. Correct normalization of the CCF requires use of the 
correct value of the mean and standard deviation for each series, 
as can be seen from the definition of the correlation coefficient r 
(Eq. fl3l])). Since AGN time series are limited in duration, as points 
near the edges drop out of the calculation, the mean and standard 
deviation of the series change; in statistical language, this means 
that the series are "non-stationary." Thus, the mean and standard 
deviation of the series needs to be recalculated for every lag, using 
only the data points that actually contribute to the calculation. In 
the case where interpolated values from one series are used, the mean 
and standard deviation should be those of the interpolated points, 
not the original points. 

(b) Significance. Once the peak value of the CCF r poa k has been found, 
we want to know whether or not the correlation found is "statisti- 
cally significant", i.e., is it likely to be real or spurious? The statis- 
tical significance of the correlation depends on the number of points 
that contribute to the calculation of r pea k, not to the total number 
of points in the series. Suppose, for example, that we have two time 
series consisting of N = 50 points, but that the maximum value of 
the CCF (say, r pca k = 0.45) occurs at a such a large lag that only 
N = 30 of the points are actually contributing at this lag. For a 
linear correlation coefficient of r = 0.45 and N = 30, the correla- 
tion is significant at about the 98% confidence level. However, if 
we erroneously use N — 50, we would conclude that the level of 
significance is about 99.9%, clearly a major difference. If we fail to 
account for the correct number of points contributing to the calcu- 
lation of r P eak and simply use the number of points in the series, we 
will overestimate the significance of the correlations we detect. 

2. Interpolation scheme. There are a couple of issues that arise in this 
regard: 

(a) Which series? In the examples shown in Fig. ^4[ we have interpo- 
lated in the emission-line light curve. Is there any particular reason 



in The Starburst-AGN Connection 2001, (Singapore: World Scientific) 44 



to choose one series over the other when doing the interpolation? In 
general there does not seem to be, unless, for example, the emission- 
line light curve is markedly smoother than the continuum light curve 
(on account of the time-smearing effect) or one light curve is much 
better sampled than the other. Usually, one computes the CCF by 
computing the CCF twicp^ interpolating once in each series, and 
then averaging the resultsEi 

(b) Interpolating function. Here we have considered only point-to-point 
linear interpolation, which results in first-derivative discontinuities 
in the light curves and the CCFs. Is there any advantage to using 
higher-order interpolation functions in the time series? Generally, 
no, higher-order functions don't improve the CCFs in any sense, and 
can be grossly misleading, as higher-order fits based on only a few 
data points become hard to control. 

3. Resolution of the CCF. Can the accuracy of a cross-correlation result be 
better than the typical sampling interval? Yes, it can, as long as the 
functions involved are reasonably smooth. The analogy that is usually 
drawn is that one can measure image centroids to far higher accuracy 
than the image size, which is true because both stellar images and in- 
strumental point-spread functions are generally smooth and symmetric. 
Statistical tests as described below suggest that uncertainties of about 
half the sampling interval are routinely obtained. 



4-5 Uncertainties in Cross- Correlation Lags 

Although cross-correlation techniques have been applied to AGN time se- 
ries for about 15 years, there is still no obvious or even universally agreed- 
upon way to assess the uncertainties in the lag measurements obtained. At 
present, the most effective technique seems to be a model-independent Monte- 
Carlo method known as FR/RSS (for "flux redistribution/random subset 
selection" )L3. 

FR/RSS is based on a computationally intensive statistical method known 
as a "bootstrap" . The bootstrap works as follows: suppose that you have a 
set of N data pairs (xi,j/j) and that linear regression yields a correlation 
coefficient r. How accurate is r? In particular, how sensitive is it to the 
influence of individual points? One can assess this by a Monte Carlo process 
where one selects at random N points from the original sample, without regard 
to whether or not any point has been selected previously. For the new sample 
of ./V points (some of which are redundant selections from the original sample, 
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while some points in the original sample are missing), the linear correlation 
coefficient is recalculated. When this is done many times, a distribution in r is 
constructed, and from this, one can assign a meaningful statistical uncertainty 
to the original experimental value of r. 

This process can also be assigned to time series, except that the time tags 
of the points have to be preserved. In effect, then, this means that redundant 
selections are overlooked; the probability that in TV selections of N points 
a point will be selected zero times is 1/e, so the new time series, selected 
at random, has typically fewer points by a factor of 1/e (hence the name 
"random subset selection"). WelsfO suggests that this should be modified in 
the sense that the weighting of each selected point should be proportional to 
i/rol, where n, is the number of times the data point (xi,yi) is selected in a 
single realization. This is closer in philosophy to the original bootstrap, but 
it has not been rigorously tested yet. 

The other part of the process, "flux redistribution," consists of changing 
the actual observed fluxes in a way that is consistent with the measured 
uncertainties. Each flux is modified by a random Gaussian deviate based 
on the quoted error for that datum (i.e., after a large number of similar 
modifications, the distribution of flux values would be a Gaussian with mean 
equal to the data value and standard deviation equal to the quoted error). 

A single sample FR/RSS realization is shown schematically in Fig. 
For each such realization, a cross-correlation is performed and the centroid 
is measured. A large number of similar, realizations will produce a "cross- 
correlation peak distribution" (CCPD)N, as shown in Fig. jH The CCPD 
can be integrated to assign formal uncertainties (usually ±ltr) to the value of 
Tcent measured from the entire data set. 

5 Observational Results 

With the background provided in the previous sections, we can discuss some 
of the more important observational results. We will begin with emission- 
line variability, since the results obtained thus far have been relatively less 
ambiguous than the results on interband continuum variations. 

5. 1 Size of the Broad-Line Region 

By the late 1980s, the potential power of reverberation mapping was widely 
understood, but the high-intensity monitoring programs required to extract 
this information had not been carried out, mainly on account of sociological 
barriers: the sampling requirements, in terms of time criticality and number 
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Figure 29. The filled circles show the Mrk 335 light curve from Fig. but now with error 
bars shown. The open circles show a single FR/RSS Monte Carlo realization; the points 
are selected from the total subset at random, and the fluxes are adjusted as Gaussian 
deviates. Note that some of the original data points are not seen because the selected flux- 
adjusted points cover them. The realization shown here gives r CC nt = 17.9 days, compared 
to Tccnt = 15.6 days for the whole data set. 



of observations, were far in excess of anything that had been done in the past 
on oversubscribed shared facilities. In late 1987, a large informal consortium 
known as the International AGN Watcf^ was formed in an attempt to obtain 
sufficient observing time with IUE and various ground-based observatories 
to measure emission-line lags for the Seyfert 1 galaxy NGC 5548, even at 
that time one of the best-studied, variable AGNs. This turned out to be 
an enormously successful projectJiSE3'E3 that provided impetus for additional 
similar projects by the AGN Watch and other groups (see PetersonEl for a 
review of the field in the wake of the first monitoring project on NGC 5548). 
Some of the light curves and CCFs from the original NGC 5548 project are 
shown in Fig. pH] . 

Since this initial successful campaign, a handful of AGNs have been well- 
monitored in the UV and optical (NGC 3783, NGC 4151, NGC 7469, 3C 390.3, 
Fairall 9, in addition to NGC 5548) and nearly three dozen Seyfert galaxies and 
low-luminosity quasars have been closely monitored in the opticalE3'E3. The 



"More information about the International AGN Watch and all papers and date 
obtained are a.va.ilah le on the AGN Watch website at http://www.astronomy.ohio- 



5tate.edu/~agnwatch / 
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Figure 30. Multiple Monte-Carlo realizations such as those in Fig. |29| are used to build up 
a cross-correlation peak distribution (CCPD). Relative to the measured CCF centroid for 
the whole data set (r = 15.6 days), the ±1<t width of this distribution is +7.2, —3.1 days. 



most important general conclusions reached from these monitoring programs 
are: 

1. The UV/optical continua vary together, to within a couple of days. 

2. The highest ionization emission lines respond most rapidly to continuum 
changes, implying that there is ionization stratification of the BLR. 

3. The BLR gas is not in pure radial motion, since no unambiguous differ- 
ences in the time scale of the response of the blue and red wings of the 
emission lines have been detected. 

Cross-correlation results show that the—size of the H/3-emitting region 



and that, in general, sizes derived from other lines usually scale as 
r(Hen, Nv) w 0.2r(H/3) and r(Lya, Civ) » 0.5r(H/3). A slightly differ- 
ent version of the BLR radius-luminosity relationship is shown in Fig. |3^, in 
which r oc i°- 6±01 . This is close to the slope of the relationship predicted 
by the most naive sort of scaling: since to lowest order AGN spectra all look 
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scales with average continuum luminosity ai 
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Figure 31. Light curves (left-hand column) and CCFs (right-hand column) for NGC 5548, 
from the first International AGN Watch monitoring project in 1988—89. From the top, the 
light curves are the UV continuum at 1350 A, the optical continuum at 5100 A, Hell A1640, 
LyaA1215, CivA1549, Cm]A1909, and H/3A4861. The continuum fluxes arc in units of 
10~ 15 crgs s _1 cm -2 A -1 , and the line fluxes are in units of 10 _13 ergs s — 1 cm -2 . The 
CCFs arc computed relative to the UV continuum at the top; the top panel is the UV 
continuum ACF. Data from Clavel et al. 14 and Peterson ct al. 67 . 
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Figure 32. The BLR radius, as measured by the H/3 lag, as a function of average optical 
continuum luminosity. The open symbols are Palomar-Green quasars from Kaspi et al. 36 , 
the filled circles are Seyfert 1 galaxies from Wandel et al. 88 . The large symbols are narrow- 
line (Vfwhm < 2000 km s" 1 ) AGNs, and the filled diamond is NGC 4051 75 . The dotted 
line is the best fit to all the data except NGC 405f . 



alike, the ionization parameters (Eq. (|T^)) and particle densities must be the 
same in all of them. Thus, inverting Eq. (p"3|), we predict that 

where in the last step we have also assumed that the shape of the ionizing 
continuum is not a function of luminosity. Despite these gross assumptions, 
the agreement with the data is not bad, although the large scatter alone tells 
us that this is not the whole story. 



5.2 AGN Black-Hole Masses 

Reverberation mapping is one method of measuring AGN central masses via 
the virial relationship 

fro- 2 

M = — > (40) 
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where / is a dimensionless factor of order unity that depends on the geometry 
and kinematics of the BLR, a is the emission-line velocity dispersion, and 
r is the size of the emitting region. Measurement of the emission-line time 
lags provides the ingredient that has been_missing since the first attempts 



estimators, advantage of using the BLR to provide an estimate of the mass of 
the central source is that it is located very close to the central source (within 
~ 10 3 i? grav ), leaving little doubt that the central mass is in fact a black 
hole. On the other hand, the kinematics of the BLR are not yet understood 
(see below), and non-gravitational forces might have a strong effect on gas 
motions. For the virial method to be applicable, the BLR kinematics must be 
dominated by gravity. Even without understanding the detailed geometry and 
kinematics of the BLR, we can test the virial hypothesis by comparing lags 
and line widths measured in a single AGN: all lines must give the same virial 
mass, even though not all the line-emitting material needs to have common 
kinematics. The three AGNs for which this can be easily tested are shown in 



Even if this were not true for all lines, it may be true for some lines, 
and a given line must always yield the same mass. Only in the case of 
NGC 5548 is there sufficient information on the long-term behavior of a single 
line (H/3) for this testjto be applied, and the data seem to be consistent with 
the virial relationshipO. We expect, then, that as the continuum brightens, 
the emission- line lag increases (see Fig. |3^), and the emission-line becomes 
narrower. This does seem to be the case. 

Virial masses based on H/3 line reverberation as a function of optical lumi- 
nosity are shown in Fig. ^5[ There is considerable scatter in the relationship, 
but it is nevertheless clear that higher-mass black holes are found in higher- 
luminosity AGNs. Some of the scatter in this relationship may be attributable 
to differences in accretion rate or radiative efficiency: the lower end of the en- 
velope, for example, seems to be dominated by narrow-line Seyfert 1 galaxies, 
which are thought to have relatively high accretion rates (and thus luminosi- 
ties) for their mass. Additional factors, such as inclination of the system, may 
also contribute to the scatter. But we are, finally, beginning to see the first 
indications of a mass- luminosity relationship for AGNs. 

5.3 Emission-Line Transfer Functions 

While existing AGN monitoring data have been of sufficient quality and quan- 
tity to obtain cross-correlation lags, recovery of transfer functions has, not sur- 
prisingly, proven to be far more difficult. Existing transfer function solutions 
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Figure 33. Line width in the rms spectrum plotted as a function of the distance from the 
central source (upper horizontal axis) as measured by the emission-line lag (lower horizontal 
axis) for various broad emission lines in NGC 7469, NGC 5548, and 3C 390.3. The dashed 
lines are best fits of each set of data to the relationship logV FWH M = a, + b log cr. The 
solid line shows the best fit to each set of data for fixed b = —1/2, yielding virial masses 
of 8.4 x 10 6 M , 5.9 x 10 7 M , and 3.2 x 10 s M Q for the three respective galaxies. From 
Peterson & Wandel 72 ©2000 AAS.. 



tend to be very noisy and ambiguous. 

In Fig. |3^, we show a sample transfer function for the H/3 emission line in 
NGC 5548; this transfer function is based on data extending over more than 
a crossing time, so the reader is cautioned against concluding too much from 
this example, since it is based on data that span a long interval compared to 
the BLR dynamical time scale (Eq. (p~9|) ) ■ The structures seen in this transfer 
function do not correspond to those seen in any of the simple models that 
we have described earlier. Note in the onc-dimcnsional transfer function the 
low amplitude of the response,_at zero lag, first noticed even with just the 
first year of AGN Watch datatll; this indicates that there is little response 
due to material along our line-of-sight to the continuum source, suggestive of 
either a low-inclination disk (i.e., there is little gas along the line of sight) or 
anisotropic line responseEiL Whether or not other lines have small response 
at small lag is less certainc3'E3; this effect may be seen clearly only in H/3 on 
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Figure 34. The H/3 emission-line lag as a function of optical continuum flux for the Seyfert 
1 galaxy NGC 5548. The solid line shows the best power-law fit to the data, r ex. F ' 4 ^' 2 , 
and the dotted line shows the best fit to the naive theoretical prediction r ac F 1 / 2 (Eq. (B9I)). 
The dashed vertical line shows the estimated starlight contribution, as in Fig. pL 



account of the large lag. 

Fig. [37] shows an attempt to recover the C IV transfer function from 39 
daily observations of NGC 5548 made wiJJaJJSTJn 1993. These data have 
been used in several independent analysesE9tJ : tirLl with no consensus on the 
interpretation, and indeed with quite contrary conclusions about the kine- 
matics: Wanders et al.Ej favor no radial motion, Done & Kroliktj find a hint 
of radial infall (also previously suggested by Crenshaw & Blackwelllia on the 
basis of the first year of IUE data), and Chiang & MurraytLd and Bottorff et 
al.tJ fit the data with different radial-outflow models. 

A crude measure of the velocity field might be obtained by cross- 
correlating parts of the emission line: for example, comparison of the response 
times for the red and blue wings could be used to detect radial inf all/outflow. 
Similarly, one could compare the response times for the line wings with the 
line core in an attempt to detect virial motion, i.e., V oc r -1 / 2 . There have 
been a number of reports that the emission lines wings respond faster than 
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Figure 35. The central mass of AGNs, as measured from ) 
of optical continuum luminosity. Symbols are as in Fig. 
AAS. 



> line reverberation, as a function 
From Peterson et al. 75 ©2000 



cordis, as expected for virial motion, but detection is always weak. 

A two-dimensional transfer function for_the Civ-Hen spectral region in 
NGC 4151 is shown Fig. [38]. Ulrich & HorneE3 argue that there is a hint that 
the red wing responds slightly more rapidly than the blue wing, which is the 
expected signature for radial infall. 

5.4 Emission-Line Profile Variability 

If the BLR velocity field is ordered (i.e., V — f(r)), then the emission-line 
profiles should vary in response to continuum changes: the propagation of 
excitation inhomogeneities through the BLR should result in irregularities in 
line profiles that evolve on time scales comparable to the emission-line lags. 
Peterson et al.Lia have examined profile variations in several Seyfert galaxies 
that were monitored fairly intensively over several years. They conclude that 
profile variations are not obviously correlated with cither the continuum or the 
total emission-line flux; profile variations are not reverberation effects. They 
also find that most H/3 profiles can be modelled fairly well with three Gaussian 
components that are fixed in location and width and vary only in relative flux, 
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Figure 36. A transfer function for H/3 in NGC 5548, based on AGN Watch data from 
1988-1996. Courtesy of S. Collier. 
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Figure 37. A transfer function for C iv A1549 in NGC 5548, HST observations 40 from 1993. 
Courtesy of S. Collier. 



which thus allows multiple-peaked profiles and variable asymmetry. They also 
conclude that profile irregularities can occur quickly, appearing on time scales 
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Figure 38. A transfer function for the C iv-Hc II region of NGC 4151 based on IUE spectra. 
The strong double-peaked appearance of the C IV line is due to a deep slightly blueshifted 
absorption feature in this line. From Ulrich & Home 85 ©1996 Blackwell Sci. 



as short as days or weeks, but they can then persist for a very long time, 
months or even years. Sample H/3 profiles for two well-studied Seyferts are 
shown in Fig. |3^. 

5.5 Multiwavelength Continuum Variability 

As noted earlier, one of the important conclusions reached from the original 
AGN Watch monitoring program was that the UV and optical continua vary 
with no apparent time delay between them, at least to the accuracy of the 
experiment ( £ 1 day in the best cases). An important consequence of this is 
that AGN continuum variability is not due either to mechanical accretion-disk 
instabilities or to variations in the accretion rate, since such effects would be 
expected to propagate through the disk (and thus across the spectrum) on 
much- longer sound-crossing (Eq._(JLl)) or drift time scales (Eq. ([i"2"|)), respec- 
tively. The suggestion was madeLLZFEj that the UV/optical variations might in 
fact be driven by X-ray variations, possibly in a manner consistent with the X- 
ray reprocessing models that were being developed to account for the 
reflection hump and the equivalent width of the 6.4 keV Fe Ka line! 
both of these features suggest that something like half of the emitted hard 
X-rays in AGN interact with "cold" (not highly ionized) matter that covers 
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Figure 39. Emission-line profiles for H/3 in NGC 5548 (left-hand column; International 
AGN Watch data) and Akn 120 (right-hand column; data from Peterson et al. 69 ). These 
are representative profiles obtained over several years of monitoring. The H/3 narrow-line 
component is particularly noticeable in the spectra of NGC 5548. 



much of the sky as seen from the source (and which might in fact be the 
accretion disk itself). If it is supposed that hard X-rays are produced above 
the disk plane near the axis of the accretion disk, then hard X-ray radiation 
striking the disk should be reprocessed into lower-energy continuum photons. 
Moreover, the UV/optical variations should follow those in the X-rays, with 
the shortest time delays for the higher-energy photons that are produced pre- 
dominantly in the central regions of the accretion disk. The temperature 
structure of a classical thin accretion disk is given by 



3GMM 
8nar 3 



R„ 



-1/2 



(41) 



in The Starburst-AGN Connection 2001, (Singapore: World Scientific) 57 



and assuming that AGNs are accreting close to the Eddington rate with about 
10% efficiency, the radial temperature structure becomes 

T{ r) - 6.3 x 10= '" (" Y 1 " (S.) K , (42) 

from which we expect a peak in the spectral energy distribution in the UV/soft 
X-rayggon, which as noted earlier might in fact be the origin of the big blue 

If we suppose that a thin accretion disk is irradiated by an X-ray source 
on the disk axis, we should see the inner, hotter part of the accretion disk 
respond before the outer, cooler parts. From Wien's Law, the radiation at 
wavelength A arises primarily at a particular temperature corresponding to a 
particular location in the disk, i.e., A cx T- 1 cx r 3 / 4 , so the difference in the 
response times at different wavelengths should be naively 

ArocA 4/3 . (43) 

It is straightforward to show that this relationship holds for an irradiated disk 
that is heated locally by hard X-rays. 

The predicted wavelength-dependent continuum time delays have been 
detected reliably in the case of NGC 7469, which was monitored intensively in 
1996Ea'E3t3. Fig. shows the result of measuring the UV and optical fluxes 
across the entire monitored UV and optical spectrum, and cross-correlating 
each of the resulting light curves with the shortest-wavelength UV contin- 
uum. Relative to the short-wavelength UV, continuum lags in other line-free 
continuum bands are detected at no less than 97% confidence throughout the 
UV/optical region. 

While the evidence shows that there are in fact wavelength-dependent 
continuum lags in NGC 7469, this is not an unambiguous detection of a clas- 
sical thin accretion-disk structure. There are several complications: 

1. Similar wavelength-dependent lags have not been reliably detected in 
othet-sources. This effect may be marginally present in the case of NGC 
4151E3, but it is not found in NGC 3516E3. However, the sampling rates 
obtained for NGC 7469 are far superior to those of most other exper- 
iments, and the upper limits on UV/optical continuum lags for other 
AGNs seem to be consistent with predictions based on (highly uncertain) 
scaling of the other sources to the NGC 7469 result, except possibly in 
the case of NGC 3516. 

2. Korista & Goad0 have suggested wavelength-dependent lags are consis- 
tent with the expectations for diffuse emission from the emission-line 
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Figure 40. Lag as a function of wavelength for NGC 7469. The average flux in each wave- 
length band is cross-correlated with the shortest-wavelength UV (1315 A) continuum, yield- 
ing a "lag spectrum". The expected dependence for an irradiated thin accretion disk, 
t oc A 4 / 3 , is also shown. The broad emission lines are prominent because they have longer 
response times than the adjacent continuum. 



clouds (Balmer and Paschen continuum). Observations over a wider 
wavelength range are needed to distinguish between an accretion-disk and 
a diffuse-continuum origin (specifically, the NGC 7469 observations did 
not cover the crucial wavelength range around the Balmer edge ~ 3650 A) . 

3. The hard X-ray flux variations in NGC 7469 are not correlated with 
those in the UV/opticalO. The hard X-ray spectrum was monitored 
with RXTE at the same time as the UV/optical monitoring campaign. 
Comparison of the light curves shows that the minima occurred at about 
the same time, but the hard X-ray maxima occurred approximately 4 
days later than the corresponding maxima in the UV continuum light 
curve. This is a prima facie argument againaLthe reprocessing scenario. 
However, a more complete spectral analysisE3 enabled by an improved 
background model for RXTE finds that the UV flux variations are cor- 
related with variations in the X-ray spectral index, in the sense that the 
X-ray spectrum is softer when the UV is brighter. This suggests that 
the soft X-rays are strongly correlated with the UV; an increase in the 
number of UV seed photons cools the Comptonizing corona above the 
accretion disk, softening the spectrum. 

By way of contrast, a recent study of high-energy variability in NGC 
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Figure 41. Comparison of hard X-ray and optical spectral variations in NGC 4051 in dif- 
ferent X-ray states. Panel (a) shows the 2-10 keV flux measured with RXTE as a function 
time for the three-year period 1996-1998. Panel (b) shows the rms optical spectrum during 
an X-ray active period in 1996, Panel (c) shows the rms optical spectrum during an X-ray 
quiescent period in 1998. Note the absence of strong He II A4686 emission in panel (c). From 
Peterson et al. 75 ©2000 A AS. 



5548 shows clear relationships between the hard X-rays and the lower-energy 
photonal3: the extreme ultraviolet (EUV) variations lead those in the soft 
X-ray region by about 3.5 hrs and those in the hard X-ray by about 10 hrs. 
The temporal order and time scales are consistent with production of both 
the soft and hard X-rays by Compton upscattering of lower-energy photons. 

If the X-rays are indeed Comptonizcd lower-energy photons, are the seed 
photons in the EUV (~ 100 A) or UV (~ 1000 A)? The only good comparison 
available is based on simultaneous UV and EUV monitoring of NGC 5548 in 
1993E3. Cross-correlation of the UV and EUV light curves shows that the UV 
leads, with r ccnt = Q-lig 2 days, which is consistent with zero lag and leaves 
the question unanswered. 

A final recent observation that should be mentioned is based on a three- 
year combined X-ray (RXTE) and optical monitoring program on the narrow- 
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line Seyfert 1 galaxy NGC 4051L3. During the first two years of this program, 
NGC 4051 behaved in a "normal" fashion for narrow-line Seyfert Is, with 
rapid, violent X-ray variability and much lower-amplitude optical variations. 
During the third year, NGC 4051 went into a very low X-ray state. In Fig. ^l], 
the top panel shows the X-ray light curve. The bottom two panels show rms 
optical spectra, which isolate the variable parts of the spectrum, obtained 
during the two periods indicated, one during a high X-ray state and one 
during the low X-ray state. They are remarkably different. The high-state 
rms spectrum (lower left) shows that the optical continuum and the H/3 and 
He II A4686 emission lines were all varying strongly. However, the low-state 
rms spectrum (lower right) shows that the optical continuum and the H/5 
emission line are still present and variable, but the Hen line has vanished 
from the rms spectrum, meaning that it is absent or constant. Whether or 
not it has completely vanished is difficult to determine on account of the 
blending of He II with strong Fe n emission in the mean spectrum, but an 
attempt to remove the Fe II emission by subtraction of a template indicates 
that most of the Hen emission must have vanished during 1998. This indicates 
that not only has the X-ray flux dramatically decreased, but the EUV flux 
(which drives the He n variations) must have also decreased by a significant 
amount. This may suggest that the entire inner accretion disk has undergone a 
transition to a low-radiation state, such as an "advection-dominated accretion 
flow" Ell. 

5.6 Conclusion 

The few examples cited in this section give an indication of the power of co- 
ordinated multiwavelength observations for understanding the physical mech- 
anisms at work in AGNs on angular scales far too small to be resolvable with 
any current or near-term technology. These are difficult programs to imple- 
ment, but the scientific return on them is potential very large, especially as 
we can build on the significant progress that has already been made. 
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